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ABSTRACT  (Ctmllnu*  an  nwn  aid*  If  naeaaaar r and  I dan  I Ity  *r  6/ocA  miAarJ 


The  Infrared  Background  Radiance  level  in  the  atmosphere  above  60  km  is 
estimated  for  various  sight  conditions  using  the  two  existing  computer 
codes,  the  FLASH  and  the  DEGGES.  The  actual  computations  were  made  by 
using  the  programs  modified  from  these  original  codes  for  improvement  of 
processing  efficiency.. 


UNCLASSIFIED 


INTRODUCTION 


Two  processes  are  responsible  for  causing  the  background  radiation 

of  the  upper  earth  atmosphere:  the  atmospheric  scattering  of  the  solar 

radiation  and  the  natural  emission  by  the  atmospheric  constituent  gas 
1 2 

molecules.  ’ Two  computer  programs  were  developed  for  theoretical 
study  of  the  atmospheric  background  radiance  level.  The  FLASH, 
which  was  developed  by  Radiation  Research  Associates,  Inc.,  Fort  Worth, 
Texas,  computes  the  level  of  the  scattered  radiation.  Its  applicability 
has  not  been  extensively  tested  in  the  infrared  region  where  the 
scattered  radiation  level  is  significantly  low.  Degges  of  Visidyne, 
Inc.,  Burlington,  Massachusetts , developed  a computer  program  which 
produces  an  estimated  level  of  the  molecular  emission  radiance  in  the 

7 

upper  atmosphere.  The  solution  produced  by  the  DEGGES  (we  refer  to 
his  program  as  the  DEGGES  in  this  report)  covers  a spectral  region  of 
1*00  a LOOO  cm  The  scattering  process  plays  a predominant  role  in 
higher  cm  ^ range,  while  tiie  molecular  emission  dominates  the  background 
radiance  level  in  lower  cm  ^ region.  Contribution  from  both  processes 
is  expected  to  equalize  in  a range  between  3000  cm  ^ and  5000  cm  ^ . The 
present  work  intends  to  make  a critical  study  of  both  programs  in  this 
region.  Tt  in  hoped  that  the  results  obtained  by  both  programs  arc 
eventually  combined  in  a proper  perspective. 

This  report  describes  the  results  which  we  obtained  during  the 
first  stage  of  our  multiyear  program.  Our  main  effort  during  this 
stage  was  to  modify  both  programs  into  a form  readily  adaptable  to  our 
computing  facility,  the  CDC  6600  operated  under  the  NOG  1 2,  in 
addition  to  familiarize  ourselves  with  their  algorithms.  Oeveral  10 
operations  in  the  original  programs  were  found  rather  inefficient  and 
time-consuming.  Significant  improvements  have  been  achieved  in  the 
overall  computation  speed  by  streuml ining  these  inefficient  operations. 
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ATMOS  PH  EH  1 0 PARAMETERS 

The  basic  physios  Involved  In  the  scattering  and  the  molecular 
emi  ion  process  are  well  Known  In  general.  Roth  programs , KnASII  and 
DLGGKS,  need  u p.'iramet  r i /ml  knowU'dge  of  the  atmosphere  In  which 
these  processes  take  place.  Some  ot-  the  atmospheric  parameters  are, 
at  present,  known  only  in  a limited  accuracy.  It  is  desirable  that 
both  programs  perform  their  computations  based  on  a unified  atmospheric 
model  by  using  an  identical  set  of  the  at mospherlc  parameters,  and 
rather  unfortunate  at.  present  that  there  exist  some  discrepancies  in 
these  parameters  used  by  them.  The  SHIRE  program  was  planned  to 
provide  us  an  improved  knowledge  of  these  parameters.  It  is  hoped 
that  the  data  analysis  to  be  conducted  ir.  this  study  would  eventually 
bring  their  improvement. 

The  temperature  profile  of  the  atmosphere  assumed  in  the  FLASH 
and  in  the  DEGGES  are  different,  as  shown  in  Fig.  1.  The  profile  may 
be  of  a secondary  significance  in  the  FLASH  computation,  while  it  is 
to  produce  a first-order  effect  on  the  computation  by  vhe  DEGGES.  No 
attempts  have  been  made  to  run  both  programs  using  a unified  profile. 

The  total  gas  density  present  in  each  layer  of  the  atmosphere  is 
nn  important  parameter  which  affects  the  accuracy  of  the  FI.ASH 
computat ion,  while  the  data  on  the  individual  constituent  species  are 
not  for  the  computation  of  the  upper  atmosphere  scattering  where  the 
molecular  absorption  is  a negligibly  small  correction  to  t he 
scattering  computation.  Fig.  shows  the  total  gas  density  assumed 
in  the  FLASH  calculation.  The  chemical  composition  of  each  species 
given  as  a function  of  the  altitude  must  be  precisely  known  in  order 
to  make  the  DEGGES  computation  realistic.  The  total  gas  density 
affects  the  DEGGES  computation  if  u significant  amount  of  the  energy 
transfers  between  the  molecules  through  collisions.  In  the  upper 
atmosphere  the  eollisional  processes  affect  little  to  the  problem  of 
the  radiative  balance.  Fig.  3 through  Fig.  tl  show  the  molecular 
density  of  various  species  assumed  by  the  DEGGES  given  as  a function 
of  the  altitude. 
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SOLAR  IRRADIANCE 

The  solar  radiation  is  the  primary  and  the  secondary  source  for 

the  upper  atmosphere  background  radiation  which  is  originated  either 

by  the  scattering  or  the  molecular  emission  process.  The  irradiance 

value  is  a parameter  affecting  the  computation  accuracy  of  both 

programs.  Fi&  9 and  Fig.  10  show  the  solar  irradiance  assumed  in  both 

2-1  2 

programs  given  as  a unit  or  watt/cm  /cm  and  of  watt/cm  /p. 
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Program  DfeXIOES 

Most  of  the  pertinent  informations  for  this  program  are  contained 
in  the  final  report  by  Degges , AFCRL~TR-71t-0606 . Even  though  the 
present  version  of  the  DPX5GES  is  somewhat  modified  from  the  original  , 
the  report  still  serves  as  a well  documented  reference.  We  shall 
mention  only  the  points  considered  important  for  our  study. 

The  entire  program  is  currently  divided  into  three  parts.  The 
first  part  (RGND)  is  the  heart  of  the  entire  computation.  The 
integrated  band  '"adiance  corresponding  to  various  viewing  conditions 
is  computed  for  the  major  atmospheric  infrared  active  constituents: 
CO^,  H.,0,  NO,  N ,0,  0H)( , and  0^.  Table  I is  a summarized  list  of  the 
vibrational  levels  considered  in  the  computation  for  these  molecules. 

The  second  and  third  parts  are  designed  to  produce  the  spectral 
outputs  as  a function  of  the  radiance  vs.  the  wavenumber.  The 
produced  results  arc  in  two  main  categories: 

(1)  exoatmosplieric  night  and  day  limb  viewing  with  various 
tangent  heights;  and 

(2)  endoatmospheri c night  and  day  viewing  with  various 
zenith  angles . 


Figs.  11  through 
spectra  shown  in 
of  ? cm”1',  which 


•tli  shew  the  results  produced  by  the  DF0GR1 . The 
these  figures  are  plotted  with  a spectral  resolution 
is  the  current  limit. 
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Table  I 

Vibrational  Levels  Included  in  the  DEGGES  Computation 


00001 

01101 

10002 

02201 

10001 

11102 

03301 

00011 

10012 

10011 


0000.00  cm" 
667.38 
1285.41 
1335.13 
1388.19 
1932.47 
2003.30 
2349.15 
3612.84 
3714.78 


0000.00  cm 

1594.74 

3151.63 

3657.05 

3755.92 

4666.72 

5234.98 


Table  I (Continued) 


n 


0200. 

or  0220 


01*00 . 
or  01*20 ) 


00000000 

00000111 

01100001 

00011001 


0000.00  cm 

588.8 


1173.0 


1758.3 

2338.5 

2223.5 
1285.0 


0000.00  cm 
1878.0 
3727.0 


000.00  cm 

101*2.1 


0000.00  cm 
1306.2 
1526.0 
3018.1* 


Kig.  13 


Day  Limb  Viewing  at  100  Km  Tangent  Height 
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Pig.  19 


Day  Horizontal  Viewing  Radiance  at  60  Km  Altitude 
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Day  Vertical  Viewing  Radiance  at  00  Km  Altitude. 
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Day  Vertical  Viewing  Radiance  at  100  Km  Altitude. 
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Fig.  P6 


Hay  Vertical  Viewing  Radiance  at  150  (On  Altitude 
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Night  Horizontal  Viewing  Radiance  at  100  Km  Altitude 
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Night  Horizontal  Viewing  Radiance  at  150  Km  Altitude 


36 


^RVE.NJMBER.f  CM*«f -1  ) ) X !0'J 


Fig.  32 


Might  Vertical  Viewing  Radiance  at  80  Kn  Altitude 
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Program  FLASH 

'ITie  program  uses  the  backward  Monte  Carlo  method  for  studying 
multiple  scattering  and  absorption  in  spherical-shell  atmosphere. 

In  region  of  the  upper  atmosphere  where  the  present  study  is  focused, 
the  background  photons  scattered  by  either  the  Rayleigh  particles  or 
the  aerosols  would  reach  the  sensor  most  likely  by  going  through 
only  a single  deflection.  A simple  calculation  can  show  that  a 
probability  of  multiple  scattering  is  quite  small  and  in  fact  the 
FLASH  confirms  this  notion.  However,  at  present,  our  calculations 
of  the  scattered  radiance  level  are  all  carried  out  by  the  FLASH. 

After  considerable  difficulty  running  the  FLASH  during  the  early 
stage  of  our  study,  even  with  the  assistance  available  from  Blttttner, 
we  eventually  have  succeeded  to  decipher  a general  structure  of  the 
FLASH.  We  believe  that  the  program  execution  devised  by  us  is 
proper.  Currently  the  program  runs  with  the  input  data  specified  by 

(1)  the  spectral  region  in  cm 

(2)  the  sensor  altitude  in  km; 

(3)  the  sensor  azimuth  angle;  and 

(U)  the  solar  angles  (a  set.  of  7 values).  (See  Appendix  A 
for  structure  of  Job  Submit  Deck.) 

The  sensor  zenith  angles  are  automatically  set  for  values  of  17  tangent 
heights  ranging  from  50  km  to  160  km  with  10  km  increment.  The  input 
data  to  the  FI.ASH  are  first  arranged  and  streamlined  by  running  a 
small  program  called  PRFLASH.  The  source  listing  of  PRFLASH  and  other 
pertinent  informations  are  also  given  in  Appendix  A. 

Figs.  35  through  37  show  the  spectral  dependence  of  the  scattered 
radiance  level  as  a function  of  the  tangent  heights.  The  figure  for 
150  km  tangent  height  indicates  the  Rayleigh  scattering  of  l/X*' 
dependency,  reflecting  the  assumed  atmospheric  model  by  the  FLASH. 

The  FLASH  computes  the  scattered  radiance  level  as  a function  of 
the  solar  angle  0,  the  sensor's  azimuth  angle  1(1  and  zenith  angle  p. 

Once  the  single  scattering  dominates  other  processes,  the  radiance 


Ul 


! 


level  is  essentially  controlled  by  the  scattering  angle  (h)  only. 

The  curves  shown  in  Figs.  38  through  Uo  are  plotted  against  the 
scattering  angle  ®,  even  though  the  original  data  are  computed 
with  various  sets  of  (6,'  $,  p).  See  Fig.  Ul  for  these  angles. 

These  curves  strongly  indicate  that  the  radiance  levels  are  indeed 
a function  of  the  scattering  angle  ® only,  and  that  other  processes 
considered  in  the  FLASH  are  negligible. 

The  computed  results  by  the  FLASH  and  by  the  DFOGES  are  super- 
imposed in  Figs.  h2  through  U5.  The  curves  by  the  FLASH  are  for  the 
scattering  angle  close  to  90°,  the  worst  expected  case  in  the 
computed  results.  The  curve  for  the  limb  viewing  at  60  km  tangent 
height  shows  that  the  scattering  and  the  emission  radiance  level  are 
about  equal.  Toward  higher  tangent  height,  the  scattering  radiation 
falls  off  more  rapidly  than  the  emission  radiance. 
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Fig.  18  Angular  dependence  of  scattering  intensity  level 

for  50  km  tangent  height  (arbitrary  intensity  scale) 
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Fig. 


39  Angular  dependence  of  scattering  intensity  level  for 
100  km  tangent  height  (arbitrary  intensity  scale) 
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Fig.  Uo 


Angular  dependence  of  scattering  intensity  level  for 
150  km  tangent  height  (arbitrary  intensity  scale) 
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Fig.  1*2 


Day  Limb  Viewing  Radiance  at  60  km  Tangent  Height 
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Fig.  Ill 


l'nv  Mini'  Viewing  Radiance  a t.  80  km  Tangent  Height 
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Fig.  t»  5 


Day  I.imh  Viewing  Radiance  at.  l*iO  km  Tangent  Height 
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1.  Source  Listing  of  PRFLASH  and  CRW 

PROGRAM  PRFLASH ( TAPE 1 * TAPES*200B » TAPE7=200B » TAPE8=200B » TAPE9 , 
1TAPE10S=200B ) 

DIMENSION  WREF ( 128 ) 

DIMENSION  A(128) ► X ( 16 ) * A A ( 16 ) 

DIMENSION  ATHETA( 16) 

101  FORMAT (8F1 0.3) 

102  FORMAT (F7. 2) 

103  FORMA  T ( F 1 0 . 4 * 1 1 0 ) 

104  FORMAT ( 15  *F8 .0  > 

103  FORMAT (8110) 

DO  25  1=1*128 
25  WREF  ( I ) -"0 . 

REWIND  7 
REWIND  8 
REWIND  5 
REWIND  9 
REWIND  1 

READ (7*101  ) W * AL  * AZ 
READ<7» 101 ) ( ATHETA ( I ) * 1 = 1 *7) 

WX  = W 

27  CONTINUE 

READ (10*104)  N * WREF ( N ) 

IF < N . NE . 70 ) GO  TO  27 
DO  28  I ==39  * 70 
L=  I 

IF ( W • LE . WREF ( L ) ) GO  TO  29 

28  CONTINUE 

29  LA=L-1 

DO  30  I =3 *35 
l.=I 

IF ( W . LE . WREF ( L ) ) GO  TO  31 

30  CONTINUE 

31  LB==L~l 
LC=0 
J=9 

DO  10  1=1*9 
X( J) =10000. /WX 
WX=WX+20. 

J=J  1 

10  CONTINUE 

XX  = 0 ♦ 5#  ( X < .1 ) +X  < 9 > ) 

WA=W  100. 

WB  WX 1100. 

CALL  CRW  < 5 > 1 ) 

CALL  CRW (2*0) 

WRITE <5* 101 ) <X< 1 ) > 1=1 *6) 

WRITE <5* 1 0 1 ) ( X ( I > * 1=7*9) 

CALL  CRW (147*1) 

WRITE (5*102) XX 
CALL  CRW (1*0) 

CALL  CRW (27*1) 

CALL  CRW (1*0) 

WRITE (3* 105)  LC * LA 
CALL  CRW (4*1) 

CALL  CRW (1*0) 

WRI TE ( 5 * 105 ) LC  * LB 
CALL  CRW (10*1) 


56 


WRITE<5> 101 ) ( ATHETA ( I ) t 1=1 ,6) 

WRITE ( 5 » 101 ) ATHET A ( 7 ) 

CALL  CRW ( 3 » 0 ) 

IC=14 

WRITE. < 5t 103 ) AL  , IC 
CALL  CRW ( 2 » 1 ) 

Y=Al.+6357 . 

Z=20 • +6357 . 

DO  20  J=1»IC 
AXX=Z/Y 

AA< J) =180.- (ASIN(AXX) *180. /3. 1415926) 
Z=Z+10. 

20  CONTINUE 

DO  21  J= 1*14 
AAA=AA( J)  . 1 
AAB=AA  < J > + ♦ 1 
AZB=AZf . 1 
AZA=AZ- . 1 

WRITE (5* 101 > AAA * AAB 
WRITE ( 5 r 101 ) AZA * AZB 

21  CONTINUE 

CALL  CRW (28*0) 

CALL  CRW (7*1) 

REWIND  9 

IF(WB.GE. 12000, ) GO  TO  60 
40  CONTINUE 

READ (9) ( A < I) *1=1*99) 

IF < A ( 1 ) . GE . WB ) GO  TO  50 

IF  < (Ad)  . GE . WA ) .AND.  ( A ( 1 ) . LT . WB ) ) 

1 WRITE  < 1 ) <A(I)*I=1*  99 ) 

GO  TO  40 
60  CONTINUE 

DO  70  1=1 r 99 
A ( I ) =0 . 

70  CONTINUE 
A < 1 ) =WA 

71  CONTINUE 

WRITE ( 1 ) < A ( I ) * 1 = 1 *99) 

A ( 1 ) =A ( 1 ) +5 . 

IF  < A ( 1 ) , LE . WB ) GO  TO  71 
CONTINUE 
50  CONTINUE 
REWIND  1 
CALL  EXIT 
STOP 
END 
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SUBROUTINE  CRW(K*L> 

DIMENSION  A ( 16 ) 

100  F ORMAT ( 8A10 ) 

DO  10  J-l »K 

READ(8> 100)(A(I)»I-1»8) 
lF(L.NE.O)  WRITE ( 5 * 100 ) ( A ( I ) » 1=1 » 8 > 
10  CONTINUE 
RETURN 
END 


Job  Submit  Deck 


JOBCARD 

ACCOUNT# 


Tl*0,  CM12000. 


ATTACH. FLASH. 

GET.FLASHDE. 

GET » FLASHB2 . 

GET . TP3DATA. 

C 0 F Y C R . T P 3 1 1 A 1 A . I A P F 1 0 . 

REWIND . TAPE  10 . 

GET  » FLASHB . 

SKIFF  . F L ASH  » 1-2 . 

COPYBF  .FLASH. TAFL9. 

PACK. TAPE9. 

COi  'YCR,  INF'U  l .TAPE  7. 

COF'YCR . FLASHHE  . TAPE 8 . 

REWIND . TAPE? . TAPES . T APE9 . 

MAP r PART. 

FLASHB2.  f OBJECT  PROGRAM  OF  PRFLASH 

COT YBF .FLASH. TAPE3. 

PACK . T APE3 . 

REWIND. TAPL1 . TAPE3. TAPES. 

FLASHB.  i OBJECT  PROGRAM  OF  FLASH 

REWIND. TAPE4. 

COF'YCR. TAFE4. OUTPUT  . 

REWIND . TAPE A . 

COF'YCR . TAFE4 » FL  ASHOT  . 
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